Statistical properties of one-dimensional random lasers 
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Statistical properties of a laser based on a one-dimensional disordered superlattice open at one 
side are studied numerically. The passive normal modes of the system are determined using the 
Feshbach projection technique. It is found that the mode competition due to the spacial hole 
burning leads to a saturation of the number of lasing modes with increasing pump rate. It is also 
responsible for nonmonotonic dependence of intensities of lasing modes as functions of pumping. 
Computed distributions of spectral spacing and intensity statistics are in qualitative agreement with 
experimental results. 
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Random lasers [1] comprise a large variety of active 
radiative systems that are based on disordered media. 
They possess one or more characteristics that distinguish 
them from the traditional lasers such as absence of the 
resonator, irregular electric field distributions and spec- 
tra, etc. A number of problems can be formulated in con- 
nection with random lasers. One group of questions is re- 
lated to light localization in disordered systems with am- 
plification and absorption [1] . Another class of problems 
deals with statistics of lasing modes [2-6], e.g., distribu- 
tion of their frequencies, intensities, etc. From a method- 
ological point of view, extension of standard laser theory 
to random lasers is an important open question, although 
some progress in this direction has been made [2, 4, 7, 8]. 

Statistical properties of random lasers were investi- 
gated in several theoretical and experimental papers. 
Experimentally, distributions of the spacial size of the 
modes and of the spacing between lasing frequencies in 
porous GaP were reported in Ref . [9] . The spacing statis- 
tics were also studied in lasers based on colloidal solu- 
tions [3], where, in addition, statistics of the emitted in- 
tensity were measured. 

Theoretical studies of this issue were mainly performed 
for lasers with chaotic resonators, which are similar, in 
certain aspects, to random lasers. In Ref. [2] random- 
matrix theory was used to describe weakly open chaotic 
cavities, where the average number of lasing modes for 
a given pump strength and the probability distribution 
of the lasing thresholds were derived. In Refs. [4-6] the 
average number of modes, its fluctuations, and spectral 
spacing statistics were computed numerically, particu- 
larly, for the case of cavities with overlapping resonances. 
Statistics of light reflection by a random laser were stud- 
ied in Ref. [10]. 

In this paper we study statistical properties of a partic- 
ular model of random lasers from first principles, without 
relying on random-matrix-type hypothesis about statis- 
tical properties of cavity modes. Regular semiclassical 
multimode laser theory [11, 12] was first applied to las- 
ing in wave-chaotic resonators in Ref. [2]. However, in the 
case of lasers based on disordered structures, this theory 



has to be modified to take into account several factors 
specific for these systems. First, in this case the open- 
ness of the system is essential and should be properly 
incorporated into the theory. In Ref. [13] the modes of 
the open cavities were considered using Feshbach projec- 
tion technique [14], and this approach was applied specif- 
ically to random laser systems in Ref. [4] . An alternative 
approach based on so-called "constant-flux modes" was 
developed in Ref. [8]. Second, the inhomogeneity of the 
background dielectric constant introduces gain-induced 
coupling between modes, which is present already in the 
linear approximation [7], but can be enhanced due to 
nonlinear hole burning effects [15]. 

In this work we rederive Maxwell-Bloch lasing equa- 
tions in the third order of nonlinearity rigorously taking 
into account the openness of the system. These equations 
are used for numerical statistical analysis of lasing from a 
one-dimensional disordered structure. We apply the Fes- 
hbach method to compute complex eigenfrequencies and 
wavefunctions of the structure under consideration. It 
is our main conclusion that the correlation of the wave- 
functions critically affect lasing properties of the system. 
In particular, they are responsible for saturation of the 
number of lasing modes and a nonmonotonic behavior of 
lasing intensities with increasing pumping. Both these 
effects are absent in the chaotic cavities described by 
random-matrix theory, in which wavefunctions are not 
correlated [2, 4]. A similar nonmonotonic behavior of in- 
tensities was recently found in Ref. [15] using a different 
approach. We also calculated statistical distributions of 
spectral spacings and mode intensities, which show qual- 
itative agreement with experiment. 

We consider a one-dimensional structure open at one 
side, which is characterized by a real non-uniform re- 
fractive index n(x) for 0<x<L, n=l for x > L, 
and an ideal mirror at a; = 0. Normal modes of this 
system have a finite lifetime even in the absence of ab- 
sorption. They can be found using the system-and-bath 
approach [13]. This approach requires division of the en- 
tire system in a closed resonator and environment, which 
interact through specified boundary conditions. In our 
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formulation wc designate the region < a; < _L as a res- 
onator with Dirichlet boundary condition at a; = L — 0, 
then the bath region x > L must be described by Neu- 
mann conditions at x = L + [16]. 

If eigenfrequencies and eigenfunctions V'a(^) 
closed resonator are known, then the true normal modes 
of the open system are obtained by diagonalizing the non- 
Hermitian matrix 

Q.{u) = 0."- mW{uj) W'^ (w). (1) 

Here fi"^ is a diagonal matrix of the eigenfrequencies lo\ 
and W(a;) is a column of coupling elements of the modes 
A with the ^-normalized bath modes 1/2/77 cos[ti;(x — L)] 
(in the system of units with the speed of light c = 1). 
Explicitly, Wa(w) = (V'A)'(-^)/v^2iwfw. The eigen- 
values of n(a;), = a;/s(a;) — iKk{co), provide fre- 
quencies and decay rates of the normal modes. Since 



n{uj) is non-Hermitian we have to differentiate between 
left, \lk{oj)), and right, |rfe(a;)), eigenvectors, which are 
biorthogonal. They can be used to define left and right 
eigenfunctions ijjl{x;uj) = J2x^^xi^)Wjk{i^)), j = hr, 
of the open system. The electric field E{x;lu) in the 
frequency representation and other relevant functions 
can be expanded in terms of these normal modes as 
E{x;oj) = Y,^Ek{uj)i)l{x]Uj). 

The gain medium is described by the polarization and 
the population difference that interact with the classi- 
cal field in the resonator. The atomic variables can be 
eliminated perturbatively from the coupled equations of 
motion [17]. Following Ref. [18] we carried out this pro- 
cedure in frequency domain, which is more convenient in 
our case than more traditional time-domain considera- 
tion. The resulting equations for mode amplitudes are 



1 f duj'dLo" 

- - U!k{ui) + iKk{u}) - ipD{uj)]Ek{u}) = -pD{uj) / D^lo') E^^iuj - uj') 

2 J (27r) ^^^^^ 

X [Bkk,k,kA^,^ - - -^") D*{-uj") Ek,{uj' - J') Eli-u;") 

+Bkkik3k2{^,^ - - w") D{-u;") Ek:„ {-i^'") E*., (-x'' - . (2) 



where the r.h.s. contains nonlinear terms in the leading 

(third) order in the field. In these equations, p is a pump- 
rate parameter (uniform pumping is assumed), D(uj) = 
[1 — i{uj — v) / ^ , Dll(a;) = (l — zcj/7||) , v is the 
atomic-transition frequency. 7^ and 7|| are polarization 
and population relaxation rates, and 

•Sfcifc2fc3fc4('^ii^2,t^3,t^4) = i / dx [V'fei(a;;wi)]* 
X Vfc2 {x; (^2) Vks (x; ^3) [Vki (x; u>4)] * (3) 

denote the overlap integrals. The electric field (in 

the time representation) is measured in units of Eq = 
■\/^7-L7||/87r (i-^i/, where d is the atomic-transition dipole 
moment. The effect of gain-induced linear mode cou- 
pling [7] was found to be small in this system and is 
neglected in Eq. (2). We do not take into account the 
effects of linear gain and nonlinearity on the lasing fre- 
quencies approximating them with solutions of equations 
uJki'jj) = ijJ, which will be denoted as ujk- In the slow- 
amplitude approximation, the mode amplitudes Ek{u) 
are assumed to be strongly peaked at uJk- The mode de- 
cay rates, eigenvectors, and wavefunctions in this case 
can be taken at the respective lasing frequencies, and 
the frequency arguments can be omitted. We transform 
Eq. (2) to the time representation and obtain rate equa- 



tions for the intensities, Ik = \Ek{t)\'^, in the form 

ik = 24 ^\Dkf - Kk -pY^Ckk'h^ , Dk = D{uJk). 

(4) 

Here terms oscillating at beat frequencies have been 
ignored. The nonlinear coupling between the modes 
Ckk' = |£)fc/|^Re [Bkkk'k'Dk] depends on the overlap in- 
tegrals (3). This is a simplified expression for Ckk' , accu- 
rate when the population inversion is time- independent. 
The correction due to population pulsations, 
ACkk' = (1/2) Re [Bkk'kk' DHcjk - tok') Dk {Dk + Dl,)] 
{k ^ k/) [6], can be disregarded for typical situation 
7|| ^ 7±- The r.h.s. of Eq. (4) has a transparent physical 
meaning of a balance between the gain, damping, and 
nonlinear saturation that prevents the intensity from 
growing indefinitely. 

The amplitude of the emitted modes outside of the res- 
onator, e™', is related to the internal mode amplitudes, 
Ek, by the input-output relation [13] with zero input: 
= —iEkit) W^rk) ■ Using this relation in com- 
bination with natural representation of matrix (1) as a 
sum of Hermitian and anti-Hermitian matrices, one can 
derive a physically transparent expression relating the 
inside and outside intensities: 

/r*=«fe4/7r. (5) 
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We model a one-dimensional disordered resonator by 
a superlattice ABAB... consisting of Ni layers of fixed 
width a. All layers A have the same refractive index ua 
and for layers B the index is drawn randomly from an 
interval (fimin, J^max)- In a-U our numerical examples the 
parameters are a ^ 1, Ni = L/a = 200, ua — 1-5, 
nmin = 0.9, and rimax = 1.3 (the refractive index is de- 
fined relative to the surrounding medium, where it is set 
to unity). This system is a periodic-on-average struc- 
ture with remnants of the bands of a periodic lattice 
with Ni/2 elementary cells. While the band gaps are 
washed out by disorder we can still define the bound- 
aries of (quasi) bands as minima of the disorder-averaged 
localization length (see inset of Fig. 2). The number of 
eigenmodes for each band fluctuates from realization 
to realization around Ni/2. We focus our simulations on 
the frequency region of the third band 2.40 < < 3.61, 
and, respectively, restrict the basis of eigenvectors V'a(^) 
used to construct matrix (1) to the modes from this band. 

We computed numerically the stationary solutions 
{Ik = 0) of Eqs. (4) which satisfy the condition Ik{p) > 
while continuously increasing pumping from zero. It 
should be understood that this method does not allow 
us to detect situations when a mode loses its stability 
without passing through zero intensity. 

The results presented below are obtained with 1700 
realizations of disorder. Figure 1 shows the pump depen- 
dence of the average number of modes (N) and its relative 
fluctuations ^var(A^)/ (N) . The atomic frequency v is 
chosen at the middle of the band and the gain width is 
equal to 20% of the band width. An important effect re- 
vealed in this figure is the saturation of (N) for large p at 
the level of 44 modes. This phenomenon is not an artifact 
of the finite basis, which contains 100 available modes. 
The saturation appears as a result of the nonmonotonic 
pump dependence of intensities Ik (p) and mode suppres- 
sion (see the inset). The nonmonotonic pump depen- 
dence was observed numerically in two-dimensional dis- 
ordered lasers as well [15]. Since in systems with uncor- 
related wavefunctions, such as two-dimensional chaotic 
resonators, where Bkkk'k' = 1 + 26kk', mode suppres- 
sion does not occur [2], we relate the origin of the sat- 
uration effect to statistical correlations between modes 
manifested via the overlap integrals Bkkk'k'- This effect 
should not be confused with saturation of the number 
of lasing modes in the regime of strong localization [19], 
where saturation is due to "exhaustion" of the number 
of available nonoverlapping localized modes. In the situ- 
ation considered here the system is in a truly multimode 
regime with overlapping modes and the saturation is a re- 
sult of a nontrivial interplay between effects of self- and 
cross-saturation. 

Before saturation the variance of the number of modes 
behaves as var(A^) oc (A^)"-^^. Deviation from this be- 
havior marks transition to the regime of strong nonlinear 
mode competition. 
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FIG. 1: (Color online) Average number of lasing modes {N) 
(solid line) and its fluctuations A/var(A) / (N) (dashed line), 
as a function of the pump rate p. For small p, {N) oc p"'^^ 
and ^var(A)/(A) oc p . Inset: intensities Jfc (p) of several 
lasing modes excited in the intermediate pump range, for one 
disorder realization. Parameters: u — 3.005 and 7x ~ 0.242. 



Figure 2 displays probability distributions of spac- 
ing Auj between frequencies of the neighboring lasing 
modes. If the atomic frequency is at the band cen- 
ter and 7_L is rather small (1% of the band width, in 
this example), then the distribution has three maxima 
(dashed line). This behavior results from the existence 
of well localized modes at the both band edges that de- 
spite experiencing a small gain still can lase due to their 
long lifetimes. As a result we have three well separated 
groups of lasing modes giving rise to three maxima in 
the spacing distribution. If the gain is centered at the 
band edge (dotted line), most of the lasing modes come 
from that edge, but modes with extremely small Kk at 
the opposite edge can be excited as well. This produces 
a second maximum at the full band width (outside of the 
plot range). 

The solid line in Fig. 2 presents local spacing dis- 
tribution for modes taken from a narrow spectral strip 
(10% of the band width) at the band edge, i/ being 
within the strip. This distribution displays mode repul- 
sion (vanishes in the limit Auj — > 0) and is similar to 
the Wigner surmise PwiAu) = {ttAuj/2) exp(-7rAw^/4) 
(dash-dotted line), which approximates the spacing dis- 
tribution in passive closed chaotic systems. The repulsion 
occurs because the length L is short enough, so that even 
states with localization length of L/10 can spatially over- 
lap. Moreover, two nonoverlapping modes may have close 
frequencies, but one of them will be localized closer to the 
opening and suppressed for moderate pump strength. 

Probability distributions for internal and output in- 
tensities of lasing modes are shown in Fig. 3. The out- 
put distribution has a singularity at zero intensity, whose 
existence is related to a very broad (several orders of 
magnitude) distribution of the decay rates of the modes 
in our system. Indeed, the modes with smaller decay 
rates nk are preferentially excited, but have lower output 
intensities according to Eq. (5). Both distributions have 
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FIG. 2: (Color online) Probability distribution of the spec- 
tral spacing between lasing modes (7x = 0.0121, p — 10~^). 
Parameters and properties for the dashed, dotted, and solid 
lines, respectively: atomic frequency f = 3.005, 2.405, 2.405; 
modes collected from the frequency interval (2.400,3.610), 
(2.400, 3.610), (2.400, 2.521); determined mean spacing ^ = 
0.27, 0.18, 0.019. Dash-dotted line: Wigner surmise (see 
the text). Inset: average localization length (in units of the 
layer width) in the third quasiband, as a function of eigenfre- 
quency lj. Horizontal line indicates current system size. 




FIG. 3: (Color online) Probability distributions P{I / {!)) of 
internal mode intensities 7^ (solid line) and output intensi- 
ties 7fc"* (dashed line). Parameters: v = 2.405, 7^ = 0.242, 

-3. 



p = 10 '^; modes are collected from the frequency inter- 
val (2.400,2.521). Mean values: {h) = 0.135, (7r*> = 
1.78 X 10~®. Insets: (a) semilogarithmic plots; (b) loga- 
rithmic plots with the fits P(x) — 6.6 x^ ''^ (internal) and 
P{x) = 0.082 ° (output), where x = 7/(7). 



approximately a power-law asymptotic behavior at large 
intensities, with different exponents (see the insets). 

It is interesting to compare our numerical results for 
the spectral spacing (Fig. 2) and the intensity distri- 
butions (Fig. 3) with available experimental data. The 
spacing distribution in porous GaP was found to be well 
approximated by the Wigner surmise [9], while in col- 
loidal solutions of Ti02 particles the mode repulsion was 
not of Wigner type [3]. At the same time, the power- 
law distribution of intensities of lasing modes found in 
Ref. [3] agrees with our simulations. 

In conclusion, we studied numerically a model of dis- 



ordered laser based on a one-dimensional open resonator. 
The passive normal modes of the system were determined 
self-consistently using the Feshbach projection technique. 
The intensities of lasing modes were found from the rate 
equations within the semiclassical third-order laser the- 
ory. Mode competition, as a consequence of the spacial 
hole burning, leads to nonmonotonic pump dependence 
of intensities and mode suppression. The number of las- 
ing modes saturates with increasing pump rate. The local 
spectral spacing distribution shows a Wigner-like mode 
repulsion. Globally, the distribution can have several 
maxima due to the quasiband structure of the spectrum. 
Distributions of the mode intensities have a power-law 
asymptotic tail. Output intensities are distributed over 
several orders of magnitude, reflecting the spread of ra- 
diative lifetimes of normal modes. 
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